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Abstract 



Birnbaum's three-parameter logistic function has become a common basis for item 
response theory modeling, especially within situations where significant guessing behavior 
is evident. This model is formed through a linear transformation of the two-parameter 
logistic function in order to facilitate a lower asymptote. This paper discusses an 
alternative three-parameter logistic model in which the asymptote parameter is a linear 
component within the logit of the function. This alternative is derived from a more 
general four-parameter model based on a transformed hyperbola. 
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Introduction 

Birnbaum (1968) introduced the three-parameter logistic (3PL) item response 
model in a contributed chapter in Lord and Novick (1968). Since then, this particular 
formulation has become a standard for investigators who wish to include a lower 
asymptote in their latent trait models. The logistic function was chosen as an alternative 
to normal ogive models (Lord, 1952), due to its more convenient mathematical properties. 

Simpler item response theory (IRT) models, such as the one-parameter logistic 
(1PL) or Rasch model (Rasch, 1960); and two-parameter logistic (2PL) model (Birnbaum, 
1957), have also been proposed, studied, and widely used. In addition, many models 
which do not conform to the common assumptions underlying most IRT have been 
investigated. Examples of these include Samejima's (1979) models for nonmonotonically 
increasing item response curves; Bartholomew's (1980) full-information factor analysis for 
multidimensional data; and Masters's (1982) partial credit model for polychotomousiy 
scored items. However, besides the three-parameter normal ogive, no other alternative 
latent trait model which is similar in shape and underlying assumptions to the 3PL has 
been extensively investigated. 

This lack of research into alternatives might be due in part to the fact that the 
3PL has served the testing community very well. The three parameters are easily 
interpretable and the model's links with the 1PL and 2PL are obvious. Also, the 
application of this model to actual test data has been facilitated by the availability of 
associated computer programs, such as BILOG (Mislevy & Bock, 1982) and LOGIST 
(Wingersky, Barton, & Lord, 1982). 
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From an analytical point of view, however, the 3PL is quite difficult to work with; 
and as a consequence, many basic properties have not to date been completely explored. 
For instance, the 3PL has yet to be shown to be (or not to be) a consistent model. While 
simulation studies have recently taken the place of analytical investigations in many 
situations, a good theoretical basis will always be preferred over a finite number of 
empirical observations. 

These analytical difficulties can also translate into practical problems. For 
example, while simultaneous confidence bands are easily found for the 1PL and 2PL given 
abilities (Hauck, 1983), this is not the case for the 3PL (Lord & Pashley, 19S8). 
Parameter estimation problems have also been experienced, especially with regard to the 
lower asymptote (Mislevy & Stocking, 1989). 

In an effort to provide an alternative avenue of research, this paper presents new 
three-parameter logistic function which possesses many of the same modeling 
characteristics of the standard (Birnbaum) 3PL but with a distinctly different 
formulation. In particular, the asymptote parameter can be written as a linear 
component within the argument of a logistic function. This alternative is derived from a 
more general four-parameter model based on a segment of a transformed hyperbolic curve. 
These two models will be referred to as the hyperbolic 3PL and hyperbolic 4PL. 
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Notation 

As mentioned above, the standard 3PL can be written as a linear transformation of 
a two-parameter logistic function. A common form of this model is 

P(8) = c + (1 - cM1.7a(6 - b)] , (1) 

where P(8) represents the probability that an examinee with ability 8 will answer a 
specific item correctly; a, 6, and c are item discrimination, difficulty, and lower asymptote 
(or pseudo-guessing) parameters, respectively; the constant 1.7 is a scaling factor; and 

Vfe) = = (2) 

1 + e z 1 + e" z 

is the logistic function. 

Associated with the logistic function is the logit transformation, denoted by \(d), 
which can be expressed as 



W) = In 



P(Q) 



(3) 



The space defined by this transformation, called the logit space, provides an alternative 
modeling environment to the usual probability correct space. 



Logit Space 

IRT models which are logistic functions can be conveniently written in the form 
P(8) = ¥1X0)]. For example, the logit transform of one form of a 2PL model is 
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W9) = 1.7a(9 - b) , 



(4) 



which is simply the argument of the corresponding ogistic function. 

Unfortunately, in the case of a standard 3PL which is a linear transformation of 
the 2PL (in the probability correct space), a logit transformation does not yield such a 
simple expression. To illustrate, the following is one form of a logit transform of a 
standard 3PL (P. W. Holland, personal communication, August 27, 1990): 



UQ) = ln{l + exp[1.7a(9 - 6) - lnc]} + lnf c ^ 

1 - c 



(5) 



This expression can be rewritten to yield an alternative formula for the standard 3PL in 
the probability correct space (C. Lewis, personal communication, August 31, 1990): 



P(9) = In 



c + exp[1.7a(6 - b)] 



(6) 



An example of a standard 3PL item response curve and its corresponding logit 
transform with associated asymptotes is shown in Figure 1. In all cases, this function 
approaches the logit transform of the 2PL (which is a straight line) in the upper tail and 
the logit transform of the c parameter in the lower tail. Also note that the curve is not 
symmetric with respect to these two asymptotes in the logit space. 



insert Figure 1 about here 
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The Hyperbolic 4PL 



There are many well-known curves within the field of analytical geometry which 
exhibit shapes similar to those of the logit transformed standard 3PL. One in particular, 
the hyperbola, also possess comparable asymptotic properties. Because of these features, 
the hyperbola was chosen as a basis for developing an alternative to the standard 3PL. 

A general equation for a hyperbola is given by 




where Z and W denote the axis coordinates; and s and r are parameters which define the 
shape of the curve. An example of a hyperbola is shown in Figure 2. 



Two transformations and a constraint are needed in order for this curve to 
resemble the logit transformed standard 3PL illustrated in Figure 1. The first 
transformation is a rotation defined by 




Z m Ycos(cx) - Xsin(a) , and 



(8) 



W - Ysin(cx) + Xcos(a) , 



where X and Y denote the new coordinate system, and 
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(9) 



Substituting (8) and (9) into (7) yields 

y = sr [y + yjX 2 I (r 2 - s 2 ) ] . (10) 



r 2 - s 2 

This manipulation rotates the curve until its upper left asymptote is horizontal, thus 
providing a lower asymptote and ensuring that the resulting model is monotonically non- 
decreasing. Note that embedded in Equation 10 is the constraint that only the upper 
segment of the curve will be considered (i.e., the part of the curve which lies above the W 
axis in Figure 2). 

The second transformation is a translation needed to center the curve properly. 
This is achieved through the following reparameterizations: 

X = 9 - h , and 

(11) 

y « we) - * . 

The result of inserting (11) into (10) is the logit of a new four-parameter logistic item 
response model: 

we) = f[(B - h) + V(e - h? + g ] + k , (12) 

f > 0, g > 0. 

Three examples of this hyperbolic 4PL, with associated asymptotes, are shown in 
Figure 3. In this case, the three curves differ only with regard to the parameter g. This 
parameter reflects how quickly the curve approaches its asymptotes. As g approaches 
zero, the curve tends toward the asymptotes and coincides with them when g equals zero. 
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Note also that in the logit space, the model is symmetric with regard to the asymptotes, 
unlike the standard 3PL. 




Relationships between the hyperbolic 4PL parameters f 9 h, and k, and parameters 
a, b y and c from the standard 3PL, obtained by equating their respective asymptotes, are 
given by 



f 



1.7 



h = 6 + 



k = In 



1.7a 



.In 



, and 



(13) 



These relationships indicate that f may be regarded as a slope parameter; h is similar to a 
difficulty parameter; and k may be thought of as a lower asymptote parameter, as it is the 
logit transform of c. As in the case of the standard 3PL, the hyperbolic 4PL approaches 
the form of the 2PL as k approaches -°o (or as c approaches 0); and is equivalent to the 
2PL when k = -<». Note that the hyperbolic 4PL parameter # is not included in the 
equations of the asymptotes, but rather, as indicated above, determines how quickly the 
curves approach their asymptotes. 
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The Hyperbolic 3PL 

In order to reduce the hyperbolic 4PL model into a three-parameter form, one 
pa-ameter, or a combination of parameters, needs to be constrained in some way. 
Consider Figure 3 once more. Note how the parameter g affects the slope of the curve in 
the probability correct space. In the case of the standard 3PL, the amount of slope is 
determined mainly by the a parameter. If the hyperbolic 3PL is to behave in a similar 
fashion to the standard 3PL, one constraint option is to let g be a function of/", and denote 
it by g(f). This is a reasonable option, as Equation Set 13 indicates, since /"may be 
regarded as a slope parameter that is related to the standard 3PL a parameter when the 
respective asymptotes have been equated. Given this constraint, the logit of the 
hyperbolic 3PL can be written as 



= Ho - h) + v/(e - kf + gin 



+ k 



(14) 



One specification which simplifies this model definition is g(f) = f' 2 . Then Equation 
14 reduces to 



(15) 



UQ) = f (0 - h) + J[f (9 - h)f + 1 + k . 
Note that this constraint, or any others considered in this paper, does not affect the 
monotonically non-decreasing nature inherent to the hyperbolic 4PL, nor does it preclude 
the 2PL and 1PL as a submodels. 

A comparison of the hyperbolic 3PL defined by Equation 15 and a standard 3PL is 
shown in Figure 4. The parameters for the hyperbolic 3PL were derived from the 
standard 3PL values using Equation Set 13, in order to ensure that both curves had the 
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same asymptotes. As was the case for the hyperbolic 4PL, the logit is still symmetric 
with respect to the asymptotes. 



: .-.Insert Fi^e 4 about hIM- 



A comparison of the two models in the more familiar probability correct space is 
shown in Figure 5, using the same parameter values as in Figure 4. While these two item 
response curves are fairly similar, a better match can be obtained by optimizing the fit of 
one to the other. This can be accomplished by sampling systematically points from the 
standard 3PL (with given parameter values) and then performing a nonlinear logistic 
regression to determine the most appropriate hyperbolic 3PL parameter values. Three 
examples of such model fitting are found in Figure 6. 

Highly discriminating standard 3PL items were found to be the hardest to 
replicate. Even in these cases, though, the two curves were found to be for the most part 
within .01 of each other. This compares favorably to the maximum difference between the 
three-parameter normal ogive and logistic functions (when using a constant of 1.7), which 
is .01. 



Insert Figure •5 .abput-'here 




p ^Insert Figure 6 about here 
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For estimation purposes, at least within a maximum likelihood framework, 
practitioners usually rely on the associated likelihood equations and information matrices. 
Those related to the form of the hyperbolic 3PL defined by Equation 15 are derived in the 
appendix. 



Other Formulations and Associated Models 

Various other hyperbolic 3PL definitions can be formulated through a variety of 
model constraints and/or reparameterizations. To illustrate this point, three different 
hyperbolic 3PL model specifications are given in this section. Note that only the logit 
formulas are given, since in all cases P(8) = ¥[^(6)]. 

Example 1. Consider the model derived from using the links to the standard 3PL 
found in Equation Set 13 in order to reparameterize Equation 15: 



W6) = ± 
2 



1.7a(8 - 6) + 



1.7a(6 - 6) - In 



+ 1 



2 



(16) 



This model formulation uses the parameters common to the standard 3PL but retains the 
shape of the hyperbolic 3PL. 

Example 2. Setting the hyperbolic 4PL parameter £ (Equation 12) to zero yields a 
connected line segment curve in the logit space, as illustrated in Figure 7. , The form of 
the model may then contain a square root or absolute value sign: 
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UQ) = f(Q - h) + s/[f(Q - h)f + k 

(17) 



= f[(Q-h) + \Q -h\] + k . 
Note that this function is not continuously differentiable at 8 = h. 



: Insert Figure 7 about here 



Example 3. The model specification given in Equation 15 can be reparameterized 
by incorporating the lower asymptote parameter k, along with the parameters f and h, 
into a new intercept term to form 



W) = p 0 + P,e + \/(p 2 + p\e) 2 + i , (18 > 

where 

p 0 = k - fh , 

(5, = f , and (19) 
P 2 = -fh . 

This parameterization yields a slightly simpler form compared to Equation 15. 

This is, of course, only a small sample of the many reformulations which are 
possible. The potential advantages and disadvantages of the various model specifications 
will depend on several factors; including estimation considerations, parameter 
interpretability, and model fit. 



Alternative 3PL 
Page 14 



Conclusions 



While there has been a recent trend toward new item types, including free- 
response and interactive formats, the multiple-choice question is still the most common 
item type in use today within standardized tests. These items typically exhibit differing 
difficulty and discrimination characteristics. Due to their constrained nature, they also 
possess a lower probability correct threshold since examinees with little or no ability still 
have some nonzero chance of responding correctly by guessing. These lower item response 
curve asymptotes may also differ across items and do not depend strictly on the number of 
choices available to the examinee. 

For test developers, investigating these three inherent characteristics of individual 
items can be very useful in the construction of tests. The standard 3PL model 
incorporates these three characteristics explicitly through its three item parameters. As 
such, the standard 3PL has become one of the most investigated and implemented IRT 
models. 

On the other hand, the particular formulation of the standard 3PL (i.e., a shifted 
2PL in the probability correct space) makes this model difficult to work with from a 
analytical standpoint. In order to provide an alternative model formulation, but with 
similar modeling characteristics as the standard 3PL, this paper has introduced the 
hyperbolic 3PL whose logit is based on a transformed hyperbolic curve. A more general 
four-parameter model was first presented, then a simplifying constraint was applied to 
produce a three-parameter version whose parameters have direct links with those of the 
standard 3PL. 
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These links between the hyperbolic 3PL and standard 3PL were established under 
conditions in which the respective asymptotes had been equated. After using these 
linking formulas, the actual curves in either the probability correct or logit space can still 
differ substantially. However, if the hyperbolic 3PL is optimally fitted to the standard 
3PL, the result is usually within a .01 difference between the function values. Often the 
two curves which lie virtually on top of each other, especially when the discrimination 
parameter is in the low to moderate range. 

In addition to being able to mimic the shapes of the standard 3PL, the hyperbolic 
3PL possess other properties in common with its predecessor. These include a 
monotonically non-decreasing nature; and having the 2PL and 1PL as submodels. 

One advantage of the hyperbolic 3PL is that its lower asymptote can be expressed 
as a linear component of the logit. As can be seen from Equation Set 28 in the appendix, 
this means that the number of examinees answering an item correctly is a sufficient 
statistic for this pseudo-guessing parameter, if all other parameters are known. 

In addition, the hyperbolic 3PL item response curve is symmetrical with respect to 
its associated asymptotes in the logit space, a feature not enjoyed by the standard 3PL. 
This particular characteristic may stabilize related estimation procedures (C. Lewis, 
personal communication, June 19, 1990), especially in the case of the lower asymptote. 
Plans for testing this hypothesis through the implementation of this model in an IRT 
computer package have been set forth. 
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Appendix 



Likelihood Equations and Information Functions 



This appendix contains the usual likelihood equations and information functions 
often used in conjunction with maximum likelihood estimation. The more general 
formulas, from which equations specific to the hyperbolic 3PL model were derived, were 
taken from Lord (1980). The form of the hyperbolic 3PL used here is from Equation 15, 
and so the probability of examinee j (with ability 8,) correctly answering item denoted by 
P;„ can be written as 



where i ranges from 1 to n items, and j ranges from 1 to m examinees. 

Case 1: known item parameters 

Consider first estimating an ability 8 from an examinee's performance on n items, 
all of which have known item parameters. Note that the subscripts has been dropped in 
this case to simplify the notation. 

The first partial derivative of P< with respect to 8 is given by 



The general form of the likelihood equation that can be solved to yield a maximum 




(20) 




(21) 



where Q. = 1 - P l% 
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likelihood estimate of 9 is 



*. - P. 



A. 



(22) 



where x tf = 1 if examinee j answered item i correctly, and zero otherwise. Using (21), the 
specific form of this likelihood equation for the hyperbolic 3PL can be expressed as 



ffc - h) 



+ l 



f(e - hf + i 



= o 



(23) 



The general form of the test information function, denoted by /„, is given by 



dP 



(24) 



M PiQi 

Again using (21), a specific form of this test information function for the hyperbolic 3PL is 



/.-E/f 



- fc,.) 2 ♦ 1 



+ 1 



PiQi 



(25) 



Case 2: known ability parameters 

Now consider estimating the parameters f, h, and k for a single item based on 
results from m examinees whose abilities are known. In this case the subscript i has 
been discarded in order to simplify the notation. 

The first partial derivatives of P } with respect to item parameters f, h, and k are 
given by 
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m - h) 



+ 1 



Vm - kf + 1 
m - *) , 3 



; and 



^ - PA 
dk ' J 



The associated likelihood equations have the general form 



E 



f \ 

X- - P 
J * J 



dp 



i = o , 



(26) 



(27) 



where £ is an arbitrary item parameter. Then the corresponding likelihood equations that 
must be solved simultaneously for f, h, and k, can be derived as 



fi Ete - % - h) 



y^~- a) 2 ♦ i 



+ 1 



= o 



■ 0 ; and 



(28) 



>1 



by using (26). 

An element of an information matrix for item parameters, denoted by 1^ has the 
following general form: 
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i aPj&j 



ft PjQj aC By 



(29) 



where £ and y are arbitrary item parameters. Again using (26), the specific elements of 
the information matrix corresponding to the parameters f, h, and k, can be derived as 



>1 



/** - E/ 2 



m - h) 



+ i 



Vm - a) 2 ♦ 1 



7 ** - E*A > 



^ - Em - k ) 



'*-£ft-*l 



Vm - a) 2 ♦ i 
/ft - ^ 



+ 1 



Vm - a) 2 



+ 1 



4* - £/ 



/"ft - h) 



Jf%- hf ♦ 1* 



+ 1 



+ 1 



PyQ, , and 



(30) 
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Figure 2. A standard 3PL curve (with parameters a = 1, b = 0, and c = .2) plotted in the 
probability correct and logit spaces. The associated asymptotes (solid lines) are included in 
the logit space. 
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Figure 2. An example of a hyperbola (dotted lines) plotted with asymptotes, which are 
dependent on the parameters s and r, within a Z and W coordinate system. 
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Figure 3. Three hyperbolic 4PL curves (with parameters f = .85, h = -.82, k = -1.39, and 
g = 1, 2, and .5) and associated asymptotes plotted in the logit and probability correct spaces. 
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Figure 4. A comparison in the logit space of a standard 3PL with a hyperbolic 3PL where 
both curves have the same asymptotes. 
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Figure 5. A comparison in the probability correct space of a standard 3PL with a hyperbolic 
3PL where both curves have the same logit space asymptotes. 
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Figure 6. Item characteristic curves and corresponding residual plots resulting from fitting 
the hyperbolic 3PL to three different standard 3PL specifications. 
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Figure 7. An example of an asymptotic hyperbolic 3PL curve (with parameters / = .85, 
h = -.82, and k = -1.39) plotted in the logit and probability correct spaces. 
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